Comment on "Percolation Thresholds in the Three-Dimensional Stick System 
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We discuss a classical mistake made in earlier publications on the stick percolation problem in 
3D for generating the right isotropic configuration of sticks. We explain the observed systematic 
deviations from the excluded volume rule. New MC simulations considered by us confirm nicely the 
applicability of the excluded volume theory. 


The three-dimensional continuum percolation problem 
of permeable sticks with the form of capped cylinders 
was considered by Monte Carlo simulations in Ref. 1. 
The authors report the dependence of the percolation 
threshold on aspect ratio and on macroscopic anisotropy, 
discussing the results from the viewpoint of the excluded 
volume theory 2 . 

Throughout the paper 1 the authors claim to obtain the 
isotropic distribution of the rods orientations by generat- 
ing their 9 and ip polar coordinates randomly with a uni- 
form distribution on the [—ir/2, tt/2) and [0, 2ir] intervals, 
respectively. Following their two-dimensional study 3 
they define the measure of the macroscopic anisotropy 
of the system as: 
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Proceeding however, in the way described above, the gen- 
erated configurations will definitely not be isotropic ones, 
although their anisotropy constant (1) will be one. It is 
easy to realize that the z axes will be a privileged one, and 
percolation in this direction reached easier than in the y 
or x direction. In order to get the right isotropic distribu- 
tion for the rods orientation, their endpoints must span 
uniformly the surface of a sphere. This can be achieved 
only by choosing the 9 angle randomly with a weighted 
distribution and not a uniform one. From the surface 
element on the unit-sphere (da = sin(9) d9 dip) it is im- 
mediate to realize that the weight-factor is governed by 
the sin(9) term. 

The mistake made by the authors does not effect the 
L « r (L the length of the cylinder and r itd radius) 
limit, considered by the authors to get confidence in their 
simulation data. However, when calculating the p c crit- 
ical density at percolation and the V ex excluded volume 
of the sticks 

V ex = (327r/3)r 3 + 8nLr 2 + AL 2 r < sm( T ) > (3) 

they calculate the average of sin(j) (7 the angle between 
two randomly positioned sticks) for the right isotropic 
case, getting < sin(j) >= n/4. Calculating < sin(j) > 
for their "isotropic" configurations the result would be 
< sin(j) >= 2 1 ix. It is even more striking that in a fol- 
lowing letter 4 , confirming also the excluded volume the- 
ory, the authors do observe the systematic deviation of 


the Monte Carlo results 1 for the isotropic case respective 
to the excluded volume rule (Fig. 2 in Ref. 4), but they 
fail in explaining it. In Ref. 4 the authors argue that the 
systematic deviation is due to the fact that much smaller 
aspects ratios are required to get the right r/L — > limit. 
The difference however, is obvious and in perfect agree- 
ment with our previous affirmations. The real value of 
< sin(^) > for the "isotropic" Monte Carlo simulations 1 
should be 2/tt, which is approximately 1.24 times smaller 
than the value used (7r/4), and the systematic deviation 4 
in Fig. 2 is just of this order in the right direction. The 
error in generating the right isotropic distribution is re- 
peated in a rapid publication 5 , where the authors study 
by Monte Carlo methods the cluster structure and con- 
ductivity of three-dimensional continuum systems. The 
simulation data for the isotropic system 1 is used in a 
series of other papers 6-10 , where some tables and com- 
parison with analytical results should be reconsidered. 

New MC simulation considered by generating the right 
isotropic configuration confirm nicely the excluded vol- 
ume theory. We will discuss our new simulation datas 
and the simulation procedure in a regular paper. 
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